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Abstract 

Motivated by recent experiments on the rod-like virus bacterio¬ 
phage fd, confined to circular and annular domains, we present 
a theoretical study of structural transitions in these geometries. 
Using the continuum theory of nematic liquid crystals, we exam¬ 
ine the competition between bulk elasticity and surface anchor¬ 
ing, mediated hy the formation of topological defects. We show 
analytically that hulk defects are unstable with respect to defects 
sitting at the boundary. In the case of an annulus, whose topol¬ 
ogy does not require the presence of topological defects, we find 
that nematic textures with boundary defects are stable compared 
to defect-free configurations when the anchoring is weak. Our 
simple approach, with no fitting parameters, suggests a possible 
symmetry breaking mechanism responsible for the formation of 
one-, two- and three-fold textures under annular confinement. 

Introduction 

Over the last thirty years rod-like viruses, optically visualiz- 
able, have been established as unique prototwe systems for ex¬ 
perimental studies of liquid crystalline order^. In recent exper¬ 
iments, /(f-viruses 0.88 /tm long and 6.6 nm thick were con¬ 
fined to wedge® rectangular® and annular® geometries. De¬ 
pending on the relative size and shape of the confining geome¬ 
try various two-dimensional (2D) equilibrium configurations of 
fd-vimses with spatially non-trivial orientational texture were 
observed. The authoi® attributes the formation of rather strik¬ 
ing annular textures with three-fold rotational symmetry to the 
finite length of fd-mds. Most theoretical work so far has em¬ 
ployed numerical Monte Carlo simulations®^, incorporating 
the range of confinement, the aspect ratio and the density of 
the constituent rods. This approach inherently depends on the 
details of the microscopic parameters. Although very useful in 
modeling specific small systems these simulations do not pro¬ 
vide a universal or complete picture of the phase diagram and 
corresponding transitions exhibited by confined systems. The 
aim of this paper is to gain a basic understanding of the symme¬ 
try selection mechanisms of nematic configurations confined to 
circular geometries. 


As suggested by experimental data®, the delicate interplay 
between boundary and bulk effects is responsible for the variety 
of equilibrium structures. Confined to thin annuli, fd-vimses 
readily satisfy the preferred tangential alignment at the inner 
and outer boundaries (planar anchoring) at the expense of bend 
elastic distortions in the bulk. For a disc, such a bend config¬ 
uration would inevitably lead to the presence of a topological 
defect of charge q = -\-\ at the center, in agreement with the 
Poincare-Hopf theorem® applied to line fields, since the Euler 
characteristic j of a disc is one. However, as analytically cal¬ 
culated below, the bulk defect is unstable with respect to two 
defects sitting at the boundary. As a result the nematic direc¬ 
tor satisfies the boundary conditions almost everywhere on the 
boundary (strong anchoring) except at a finite number of points, 
similar to the studies of Langmuir monolayers®^ In case of a 
finite/weak anchoring, the preferred alignment at the boundary 
can be violated over extended region. Here we show that even 
for an annulus boundary defects become energetically favored, 
though no defects are required by topology (X = 0). By in¬ 
serting pairs of positive (negative) charges at the outer (inner) 
boundary of an annulus we screen the curvature of circum¬ 
ference and ‘unbend’ otherwise bend nematic texture, similar 
to the stabilizing effe ct of the Gaussian curvature in the two- 
dimensional manifoldlSEni. Alternatively, introducing boundary 
defects into the uniform state (no bend or splay elastic defor¬ 
mations), we can curve the director field to satisfy the boundary 
conditions and lower the anchoring energy. Therefore, when¬ 
ever bulk and boundary energy contributions become equally 
important, we expect the equilibrium 2D textures to encompass 
boundary defects. 

Interestingly, the phenomenon of the transformation between 
surface and bulk defects was studied in 3D nematic droplets^. 
When the boundary conditions, set by the temperature, change 
from the homeotropic to planar anchoring, the nematic direc¬ 
tor deforms from the ‘hedgehog’ configuration (point defect 
inside) to the pair of ‘boojums’, sitting at the surface. In experi¬ 
ments with the /c/-virus®, the nature of the boundaries does not 
change, while it is the size and the topology (disc or annulus) of 
the 2D confinement which may influence the distribution of the 
topological charge at the boundary. To quantify the effect of the 
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confinement and weak anchoring conditions, we adopt the con¬ 
tinuum theory of nematic liquid crystalsE^l, assuming the one 
elastic constant approximation, such that bend and splay elastic 
constants are equal K as suggested byElfor /t/-virus. We quan¬ 
tify the relative stability and the equilibrium number and charge 
of defects at the boundary, which is governed by the ratio of the 
anchoring extrapolation length La=K/Wa (Wa is the anchoring 
strength) and the system size R. With no htting parameters we 
capture some features of experimentsSl on /^f-viruses, namely 
the appearance of one- and two-fold textures, in particular. Our 
analytic calculations for different geometries suggest the value 
of the anchoring extrapolation length La— 5 fim. 

In the following we hrst formulate the problem and charac¬ 
terize topological defects in the bulk and at the boundary. Next 
we consider special solutions for nematic conhgurations, con- 
hned to the disc and to the annulus, allowing for the presence 
of defects at the boundary. This approach is merely an ansatz, 
the limiting case of a general minimization problem. 

Formulation of the problem 

Nematic liquid crystals are usually described by the unit vector 
n, known as the director, with n = — n to obtain an orientation 
rather than a direction. In the continuum description the elas¬ 
tic Frank free energylSlis quadratic in gradients of n, including 
bend V x n and splay V • n terms in 2D. In the one elastic con¬ 
stant approximation (equal bend and splay contributions) the 
elastic energy simplihes toCZI 

^ei = ^ JJ^dxdy\Vn\^, ( 1 ) 

where the integration is performed over the domain H. In an in- 
hnite system the ground state corresponds to the uniform direc¬ 
tor held, n = const. When conhned, liquid crystals reorient to 
satisfy the boundary conditions, resulting in spatial variation of 
the director held, parametrized in Cartesian (x,y) coordinates 
by (see Fig.[T^) 

n = cos 0(x,y)ejc + sin0(x,y)ey. (2) 

Minimizing the free energy ([T]) leads to the Euler-Lagrange 
equation 

dxxd + dyyd =0, d{z)= 52lm{ log(z - Zi)‘^' }, (3) 

i 

where z = x + iy. The equilibrium conhgurations described 
in Eq. (|^ correspond to a set of defects of strength (topolog¬ 
ical charge) qt located at positions z,. Note that in the far 
held (|z —Zi| —oo) the director is uniform, which is a good 
approximation for systems with inhnite size. The total topo¬ 
logical charge in ([^, associated with the line held n in 



Figure 1 Schematic illustration of (a) a uniform line held n = 
conhned to a simple region Q with X = I ([^, where we have chosen 
as positive the counter-clockwise orientation of the boundary dQ, 
with the unit normal V pointing outside fi; (b) two kinds of 
topological defects: bulk defects (circles of charge q, (|^) and 
boundary defects (half-circles of charge mj = Adj/{2n)), related via 
the charge conservation law, Eq. 0. 

2D, is dehned by enclosing the defects at positions Zi by an ar¬ 
bitrary cont our Y\ (see Eig. [^) and computing the following 

integral®^ 

j ds{nx dsn) = j dO^J^qi. ( 4 ) 

The effects of conhnement may be studied by adding an 
effective surface energy favoring planar anchoring in the 
Rapini-Papoular form^I^I 

= ^ ( ds(n-v)'^, n-v = cos(0-^) (5) 

2 Jda 

where Wq > 0 is the anchoring strength, s is a curvilinear co¬ 
ordinate of the boundary dLl and V = cos^(s)e,: + sin^(s)ey 
(see Eig. is a unit normal to the boundary. The anchoring 
length La = KjWa is the length scale over which the director 
reorients to align tangentially with the boundary. Any simply 
connected domain fi with boundary is homeomorphic to a 
disc, whose Euler-Poincare characteristic % = 1® with 

ds{vxdsV) = (f d(p. (6) 

27: Jdn 2 k Jdo. 

We have chosen the counter-clockwise orientation of dQ. as a 
positive one (Pig. [^l. The Euler-Poincare characteristic of an 
annulus is zero (x — 0). This can be shown by i) integrating 
over two circles with opposite orientation, connected by a cut or 
ii) triangulating an annulus and counting the number of vertices 
V, edges E and faces F, yielding x — V — E +F — 0. 

If the characteristic size of our system R is much smaller 
than the anchoring length, R < La, the director remains uni¬ 
form, as shown in Pig.[^. Very strong anchoring Wa (R ^ La), 
on the contrary, forces a preferred orientation of n relative to 
the normal V (in the case of planar anchoring n _L V). Topolog¬ 
ical defects of charge qi then appear in the bulk, with qt = X 
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(Poincare-Hopf theorerrJSl). On the other hand, the most com¬ 
mon situation in experiments corresponds to R > La. In this 
case surface anchoring and the energetic cost of bulk direc¬ 
tor deformations compete with each other and bulk defects can 
be pushed towards the boundary, resulting in director deforma¬ 
tions that satisfy the boundary conditions almost everywhere on 
do., except at a hnite number of points The repulsive na¬ 
ture of the pairwise defect-defect inte racti ons favours maximal 
separation between same-sign defects^^El while anchoring me¬ 
diated interactions between topological defects have not been 
analyzed in detail before. 

To characterize defects localized on a ID boundary dD., 
we introduce their topological charge nij related to the angle 
dehcit of the director field n, which rotates in the same (op¬ 
posite) sense as the interior contour ji (see Fig. [^), yielding 
a positive (negative) charge mj — A9j/(27z). To establish the 
connection between the strengths qi and nij of bulk and surface 
defects, the topology of the system as described by j, and the 
anchoring n • V we derive th e cha rge conservation law, follow¬ 
ing the approach outlined indinifor 3D nematic systems. We 
choose a closed curve 72 going around the boundary defects mj 
at positions (see Fig. j^). Since there are no special points 
of the director held n between the curve 71, enclosing defects 
in the bulk, and ')^, we get 

27r ^ q'; — (J) ds (n X d^n) = J) ds {n x dsn). (7) 

,■ hi Jyi 

Decomposing the integral over the closed contour yi into the 
contribution for the portion around the boundary defects, given 
by — Y,j nij, and the integral along piecewise regular curves, we 
hnd 

/ ds\{nla-\-n\)dsm-\-ns,dsnT: — nT:dsnhv XX. (8) 

'--- 

1 dsQ — ds^ 

The integral of the hrst term is 2nx The other terms 
can be computed using the parametrization n — rtyV + 11 ^ X in 
the local system of coordinates and the relations dgV = dstpx, 
dsX = —dstpv. Replacing ([^ in Q, we get the following charge 
conservation law 

In Fig. we show the transformation of a conhguration of the 
director held n containing bulk defects with qi = to 

conhgurations with boundary defects with nij = + j, + j. By 
choosing an arbitrary curve 72 enclosing the defects, one can 
show that (|^ holds in all the cases considered in Fig. In fact 
the bulk and boundary defects are in the same topological class 
since the latter are obtained simply by pushing a bulk defect to 
the boundary. Here we restrict ourselves to smooth boundaries 
and thus m = ^/2. 


To hnd nematic conhgurations minimizing the total free en¬ 
ergy one needs to solve the Euler-Lagrange equa¬ 

tion ([^ in the domain fi with the boundary condition, arising 
from the vanishing of the hrst variation of + .^a) = 0 on 
dLl, given by 

/Tv-V0 —W'asin(0 —^)cos(0 —=0. (10) 

Instead of solving this variational problem numerically, we pro¬ 
pose a plausible ansatz for the angle 9 that minimizes 
We then seek approximate solutions compatible with the one-, 
two- and three-fold symmetries observed in experiments^^! and 
compare their relative energies to determine the most favorable 
conhguration. This allows us to treat the problem analytically. 
We consider conhgurations with different number of defects 
(up to 6) and topological charge \m\ ^ 5 (or \q\ ^ 1), conhned 
to a disc and an annular geometry. Since the boundaries of both 
disc and annulus have constant curvature, the position of de¬ 
fects correspond to the furthest separation along the boundary. 
The defect charge is, however, not known a priori. Using this 
approach, we compare the energetics of nematic liquid crystals 
conhned to a disc or annulus and study the interplay between 
the anchoring extrapolation length La, the system size R and the 
core size of defects e in controlling the lowest energy conhgu¬ 
rations. 

Nematic confined to a disc 

The vector held n (|^ shown in Fig. |^,e without conhnement 
can be written explicitly as 

0(+2)(2)=Im{log(z2-/;2)}^ z=re''P, (11) 

0(+4)(^) _ ^ _l_Im{iog\/P~^}. (12) 

It accounts for the pairs of topological defects at positions 
Zi = ±R with charges ^, = -1-1 and qi = and satishes the 
Euler-Lagrange equation Q. Conhning the director helds 
given by Eqs. ( [TT] i, ( [T2I 1 to a disc of radius R (Eig. [^,f) leads 
to boundary defects characterized by m = and m = re¬ 
spectively, while the director n remains almost uniform in the 
bulk. According to Eq. (|^ we expect no anchoring contribution 
for a pair of m — where the boundary condition is not satis- 
hed just at the defect core and the total charge equals j, while 
for a pair with m — +^ the deviation from the preferred anchor¬ 
ing orientation extends over a hnite portion of the boundary, 
yielding a non-zero anchoring contribution. 

Below we quantify the energetics of the nematic conhgura¬ 
tions (Eig.[^,f) and analyze the stability of conhgurations with 
boundary defects relative to the defect free, uniform state. We 
assume that the size e of the defect core satishes e ^ R and 
that the associated core energy is much smaller than the elastic 
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Figure 2 (a)-(c) Shows the transformation of a hulk defect of charge 
+2 into two +1 defects, which subsequently approach the boundary 
of the circular region, where they can be classified as charge m = + j 
boundary defects, (d)-(f) Similarly, a charge g; = +1 defect is split 
into two qi = +j defects which in turn yield two charge m = + ^ 
boundary defects. This process illustrates the construction of our 
ansatz. The parametrization of the director field configuration of (c) 
and (f) is given by o and 0, respectively. 


energy in the bulk and thus can be neglected. First we compute 
the total free energy of configuration with m = + ^. Substitut¬ 
ing the nematic angle given in Eq. ( [T^ into the an¬ 

choring and elastic free energy and ^ei given in Eqs. 0 
and Q and integrating over the polar coordinate (p, we find 


^ ^ (13) 
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where we have truncated the converging power series in e/R. 
Including the next order corrections may be appropriate for sys¬ 
tems with relatively large defect core size e. In conventional 
nematics e is of the order of nanometers, which is the charac¬ 
teristic scale of the constituent molecules. In Ref.® the authors 
assumed e = 0.88 fim for/d-virus and estimated the anchoring 
extrapolation length La ~ 1.4 /dm. 

Nematic configurations with a pair of boundary defects of 
charge m = -l-j, in a disc (Eig.|^) have zero anchoring energy, 
except for a small region 0{e) around the defects where the an¬ 
choring condition is not satisfied. Their elastic energy is, how- 

ever, four times higher than ^ ([T^ since it is proportional 
to the square of the charge. Notice also that the structure with a 
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Figure 3 Phase diagram of nematic film confined to a disc of radius R 
showing the regions of parameters corresponding to the energy 
minimum for three configurations: a uniform state (green) and the 
two defective configurations with charges m = + j (blue) and 
m = + ^ (pink) shown in Fig.|^,f. The horizontal axis is the 
dimensionless parameter R/La = RWa/K; the vertical axis is the ratio 
of the radius R to the defect core e Cl- The circles show the location 
of the crossover between the m= ^ defect configuration and the 
uniform state for La = 1 /im and various values of the radius of the 
disc (R = 5,10,15... /im from left to right). The squares denote the 
corresponding crossover between the two defective configurations for 
La = 5 /tm In both cases we have chosen e = 0.5 /im. 


bulk defect qi = +l (Fig.[^) confined to a disc is unstable with 
respect to configuration (c), since 2;rL'log2, ne¬ 

glecting the contribution from the defect core. The fact that no 
bulk defects were observed in experiments with /tZ-virusI^ sug¬ 
gests that these systems are characterized by a finite anchoring 
strength Wa- Finally, pushing defects of hi^er charge, such as 
q=+2 (Fig .|^), to the boundary of a disJ^ gives higher elas¬ 
tic energy ~ 27r/riog(2R/e) compared to the configurations 
shown in Figs. I3c,d). 

To quantify the role of anchoring L„ = K/Wa and system 
size R we compute the phase diagram. In Fig. we compare 
three nematic states confined to a disc of size R: a uniform state 
n = const, and two nematic configurations with a pair of bound¬ 
ary defects separated by 2R and charge m = -f j and m = + ^, 
respectively (see Fig. |^,f), ( fT^ . The relative energy of these 
configurations is controlled by the interplay of the defect core 
size e, the anchoring strength La and the radius R of the disc. 
Note that the core energy for a pair of defects is of the order of 
nKq^^^, independently of e, and accounting for this contribu¬ 
tion will lead to a slight shift of the coexistence curves in the 
phase diagram without changing the basic picture. The uniform 
state is stable for small systems {R) and weak anchoring (large 
La > R). A small value of e increases the region of stability of 
the uniform state. For stronger anchoring (or smaller system 
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size) the lowest energy state is one with a pair of boundary de¬ 
fects. The m — +j defects are favored for large values of the 
core size e while the configuration with m = +^ is prefetTed for 
smaller core sizes. The main result is that by tuning the system- 
size R one can drive a transition between a uniform state and a 
defective configuration selected by the minimization of the an¬ 
choring and elastic energies. The transition is controlled by two 
dimensionless length scales R/e and R/La, where e and La are 
determined by intrinsic physical and chemical properties of the 
system. 

Next, we compare our analytic predictions with experimen¬ 
tal data for the /^f-virus® confined to a disc. For weak an¬ 
choring {La = 5 fim) we find that increasing the radius R = 
5,10,15... /tm for a chosen core size e = 0.5 jim yields a tran¬ 
sition from a uniform to a ‘defective’ state with m = ^ (squares 
in Fig. [^. For strong anchoring {La — 1 jim) the transition 
is between the two defective configurations (circles in Fig. [^. 
In both cases the critical system size where the transition oc¬ 
curs is in the range ~ 5 — 10 jUm, which is compatible with 
experiments®. On the other hand, since the defect charge asso¬ 
ciated with reorientation of the fd virus was not extracted from 
the experimental data®, we cannot determine the value of the 
anchoring extrapolation length La based solely on the results 
shown in Fig. 

In the following we adapt simplified analytical model to a 
viral nematics confined to an annular geometry. Our goal is to 
gain insight into the mechanism responsible for the selection of 
the symmetry of the boundary defect arrangement (one-, two- 
or three-fold, as seen in experiments®), as well as their topo¬ 
logical charge, and estimate the value of the anchoring strength 
La within our approximations. 

Nematic in annular geometry 

Now we consider a nematic liquid crystal confined to an annu¬ 
lus of inner radius Rt and outer radius R. This geometry does 
not require the presence of topological defects since the Euler 
characteristic of the annulus is J = 0. In this case the defect- 
free ground state shown in the top inset of Fig. satisfies the 
tangential boundary condition. We will refer to this structure 
as the bend configuration. In Fig. we compare the total free 
energy ,^tot = of a uniform state and a bend state as 

a function of R/La and ^ = Ri/R G (0 : 0.7]. For small radius R 
or weak anchoring (large La) a uniform state with zero elastic 
energy and anchoring energy ^a = -|-^)/2 is ener¬ 

getically favored, consistent with the phase diagram in Fig. 
For large a bend state with no anchoring contribution and 
elastic energy = —Kn\og^ becomes energetically favor¬ 
able. According to experiments with/t/- virus®, the transition 
between these two states occurs at ~ 5 fim and ^ ~ 0.3 — 0.5. 
Using this value, we estimate a value La ~ 2.7 — 5.5 fim, which 



Figure 4 The coexistence curve of the hend state with 
n = — sin 9 Cx -h cos (p and the uniform state with n = in an 
annulus of inner radius Ri and outer radius R (see Fig.[T]for the 
definition of the coordinate system). The horizontal axis is the system 
size R scaled hy the anchoring extrapolation length La = K/Wa, the 
vertical axis is the ratio ^ = Ri/R. The hend state {,^gi > 0 0 ) is 
energetically favored above the curve, while the uniform state 
(■^a>0 0) has lower free energy in the region helow the curve. 


is of the same order as La ~ 1.4 /im obtained for/ii-vrrus in 
rectangular geometries® with the same physical properties of 
the boundary. 

We now examine a number of nematic configurations with 
k-fold rotational symmetry, with k = 1,2,3. These can be ob¬ 
tained by considering different numbers of defect pairs of posi¬ 
tive and negative charge sitting at the inner and outer bound¬ 
aries of the annulus. The total charge is conserved since 
Z = O(0. We consider the textures shown in Fig.0and evalu¬ 
ate the corresponding energies. Configurations |1] to [3] can be 
obtained by starting from the bend state and introducing nega¬ 
tive charges m = —^ at the inner boundary and positive charges 


Table 1 The total anchoring, ,^a, and elastic, contrihutions to 
the free energy (integrated over the whole system) calculated to 
lowest order in e/R for the configurations |T|-[^ (Fig.0. The director 
n is parametrized hy the angle 6 = arg [ , where jS = 1,2,3 

for the first, second and third row, respectively, and ^ = Ri/R. 
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Figure 5 Nematic textures confined to annular geometry. The size of 
dots describes the magnitude of the defect charge: m = | ([iHl] ®nd 
p7]-[^), m = 2 ([5-[§); the color corresponds to the sign, with red 
denoting positive defects {m > 0) and blue negative ones {m < 0). 

The textures have the same anchoring energy but higher elastic 
energy than |T|-[3l, hence are not included in the stability phase 
diagram shown in Fig.j^ 


m = +1 at the outer boundary of the annulus and cost less elas¬ 
tic energy from director distortion than \2] to [^, but have a high 
cost in anchoring energy compared to 0 to Note that tex¬ 
tures [2 to [U and to have the same anchoring energy (see 
Table lb. Conversely, the configurations with \m\ = j boundary 
defects ([4]-[^ cost more elastic energy associated with curva¬ 
ture of the director field than |Tl-[3], but have lower anchoring 
energy, as shown in Fig. |^. Charge conservation as given in 
Eq. ([^ requires that for all textures considered the anchoring 
energy density is the same at the inner and the outer boundaries 
of the annulus. For the configurations [4]-^ the anchoring en¬ 
ergy is simply Waf?(l -|-^)7r^^"/2, withn = 1,2,3, respectively, 
for [4]-[^ representing the number of boundary defect pairs. 
In the following we exclude the configurations |4]-[^ from our 
analysis, because their elastic/total energy is much larger than 
that of as well as of the defect free states. 

In Fig. 1^ we compare the energies of the various nematic 
textures in the annulus. In textures 23^1] the m = bound¬ 
ary defects have high anchoring energy in fat annuli (small ^) 
but lower elastic energy than the bend state over a wide range 
of parameters ^ and R/e, as shown in Fig. |§> which displays 
the region of stability of the various textures. Thus, for small 


Figure 6 Energetics of the six configurations C0-® shown in Fig.js] 

(a) The normalized anchoring energy in units of R/La as a function of 
= Rj/R. The corresponding expressions are given in Table[^ 

(b) and (c) display the region of parameters where various nematic 

textures minimize the free energy (b) is obtained by comparing the 
elastic free energy of textures 2]-^ and of the bend state (not that for 
these configurations the anchoring contribution to the free energy is 
zero or negligible); (c) is obtained by comparing the total free energy 
of textures and of the defect free configuration shown in Fig.W 

for La = 1 pm and e = 0.5 pm. 


system size, or weak anchoring, where the elastic energy dom¬ 
inates, we expect the equilibrium textures with one-fold sym¬ 
metry (|T|) to be more favorable in thin annuli and textures with 
two-fold symmetry ([^ to be favored in thick annuli. The con¬ 
figuration [ 3 ] may also be energetically accessible in a narrow 
range of parameters, consistent with the fact that it is rarely 
observed in experimentsEl. In all cases boundary-stabilized de¬ 
fect textures are the ground states of confined nematics when 
R/La^O{l).lnFig.^, assuming weak anchoring La = 7 pm 
and e = 0.5 pm, we illustrate with different colors the minimiz- 
ers of the total free energy + ^a- Defect sttuctures with 
m = are energetically favored, and therefore support our 
previous estimate of the anchoring extrapolation length. Note 
that ground states (energy minimizers) with crystalline order on 
the surface of an embedded torus are also characterized by the 
presence of positive disclinations on the exterior of the torus 
and negative disclinations in the interioi®!. 
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Concluding remarks 


References 


We have shown that, in contrast to what speculated in earlier 
worksl2H3 a number of defective textures observed experimen¬ 
tally in nematic films confined to circular and annular geome¬ 
tries can be accounted for within continuum liquid crystal the¬ 
ory. By examining the energetics of textures with localized de¬ 
fects at the boundary, and exploiting conservation of topologi¬ 
cal charge, our work captures the main features of experiments 
in /(f-virus^land provides an estimate for the anchoring extrap¬ 
olation length in these systems - L„ ~ 5 /im. 

Several challenging questions remain unanswered. Our re¬ 
sults suggest that the selection of the symmetry of observed ne¬ 
matic textures is controlled a single dimensionless parameter, 
WaR/K, and the conservation of topological charge rather than 
by details of microscopic interactions. But what is the range of 
validity of the continuum theory, i.e., is it valid when the size 
of confinement approaches a few molecular sizes? Another im¬ 
portant question concerns the use of the one-elastic constant 
approximation, which is reasonable for semiflexible polymers, 
including the fd-vims with its persistence length of the order 
of the polymer length^EI but is not accurate for rigid rods^^. A 
significant elastic anisotropy was measured experimentally for 
systems composed of tobacco mosaic virus'^ with bend con¬ 
stant larger than splay constant, \1K\. On the contrary, 
according to theoretical predictionsl^the splay elastic constant 
K\ diverges as the length of the molecule while the bend elastic 
constant depends on the rigidity of the molecules. Thus for 
(semi)flexible polymers one expects Ki > Kj. It would there¬ 
fore be of interest to quantify the energetics of topological de¬ 
fects in case of strong anisotropy of the elastic constants. The 
authorsC^ analyzed the influence of splay-bend anisotropy on 
the formation of fiber texture in discotic liquid crystals with 
fixed boundary conditions and topological defects in the bulk. 
Our approach could also be generalized to other planar geome¬ 
tries with non-monotonic curvature of the boundary, such as 
the square plates examined in^El. Finally it is well known that 
in three dimensions topological defects affect the shape of ne¬ 
matic droplets^^HI It would similarly be interesting to examine 
the interplay between defect textures and shape in thin nematic 
Aims confined by a deformable boundary. 
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